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element  mesh  and  efficiently  analyze  the  structure.  Techniques  to  automatically  generate  a  detailed  model  of  the 
structure  and  perform  required  analysis  and  design  calculations  for  the  individual  concrete  slabs  under  all  applied 
loads  are  illustrated,  along  with  methods  for  identifying  and  reporting  details  of  the  design  for  each  slab  under  the 
worst-case  loading. 

The  two-way  slab  design  approach  (based  on  American  Concrete  Institute  guidance,  ACI  318-02)  is  summa¬ 
rized,  along  with  a  description  of  methods  to  calculate  shear,  moment,  and  thrust  throughout  the  shell  elements  and 
at  selected  cross  sections  of  the  superelements.  The  design  section  includes  the  approach  used  to  modify  design  pa¬ 
rameters  based  on  analysis  results.  Finally,  guidelines  are  provided  for  modeling  the  response  of  soil  and  piles  used 
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Conversion  Factors,  Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

meters 

inches 

25.4 

millimeters 

kips  (force) 

4.448222 

kilonewtons 

pounds  (force)  per  square  foot 

47.88026 

pascals 

pounds  (force)  per  square  inch 

6.894757 

kilopascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  meter 

square  feet 

0.09290304 

square  meters 

square  inches 

6.4516 

square  centimeters 
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1  Introduction 


This  report  is  a  theoretical  manual  that  describes  the  development  of  “design 
by  analysis”  methods  for  innovative  navigation  structures.  The  modeling 
approach  is  designed  to  facilitate  the  development  of  a  complex  three- 
dimensional  finite  element  mesh  and  efficiently  analyze  the  structure.  Techniques 
to  automatically  generate  a  detailed  model  of  the  structure  and  perform  required 
analysis  and  design  calculations  for  the  individual  concrete  slabs  under  all 
applied  loads  are  illustrated,  along  with  methods  for  identifying  and  reporting 
details  of  the  design  for  each  slab  under  the  worst-case  loading. 

This  manual  begins  with  a  description  of  assumptions  made  concerning  in- 
the-wet  construction  techniques  and  how  these  assumptions  are  accounted  for  in 
this  development.  The  section  on  modeling  considerations  documents  the 
methods  used  to  replicate  typical  details  of  the  float-in  segments.  The  finite 
element  formulation  for  the  two  element  types  is  then  presented.  Shell  elements 
are  used  to  model  the  individual  concrete  slabs.  These  are  connected  with 
superelements  derived  to  provide  an  accurate  solid  model  of  the  joint  details, 
including  tapers  between  slabs.  An  innovative  approach  to  efficiently  calculate 
the  stiffness  matrix  and  the  slave-master  transformation  matrix  for  an  arbitrary 
joint  geometry  is  introduced. 

The  two-way  slab  design  approach  (based  on  American  Concrete  Institute 
requirements,  ACI 318-02)  is  summarized,  along  with  a  description  of  methods 
to  calculate  shear,  moment,  and  thrust  throughout  the  shell  elements  and  at 
selected  cross  sections  of  the  superelements.  The  design  section  includes  the 
approach  used  to  modify  design  parameters  based  on  analysis  results.  Finally, 
guidelines  are  provided  for  modeling  the  response  of  soil  and  piles  used  in  the 
foundation.  Design  algorithms  are  presented  in  Appendix  A. 
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2  In-the-Wet  Construction 
Approach 


Traditional  methods  to  construct  massive  concrete  navigation  structures  have 
required  the  placement  of  a  cofferdam  to  dewater  the  construction  site.  Several 
Corps  of  Engineers  projects  are  currently  in  various  stages  of  design  and 
construction  using  new  methods  in  which  a  concrete  shell  is  constructed  in  a  dry 
dock,  floated  to  the  placement  site,  and  lowered  on  a  prepared  underwater 
foundation.  This  “in-the-wet”  construction  approach  will  save  much  of  the  time 
and  expense  of  cofferdam  construction. 

The  fabrication  of  a  concrete  shell  that  is  light  enough  to  float  and  strong 
enough  to  resist  substantial  pressure  from  water  and  tremie  concrete  during 
construction  requires  careful  planning  and  design.  The  development  of  “design 
by  analysis”  (DBA)  methods  is  based  on  the  approach  being  used  in  the 
Pittsburgh  District’s  Braddock  Dam,  which  reflects  experience  gained  in  other 
applications  by  the  consultants  on  the  project.  The  following  design  assumptions 
are  used  in  the  DBA  application. 


2.1  Shell  Layout  Assumptions 

The  floating  structure  is  primarily  reinforced  concrete.  A  local  coordinate 
system  is  defined  so  that  concrete  diaphragm  walls  are  arranged  along  section 
lines  in  both  the  E-W  and  N-S  directions.  The  local  coordinate  system  is  defined 
such  that  water  flows  through  the  structure  from  north  to  south.  Even  though  the 
height  and  thickness  of  the  walls  may  vary  along  the  length  of  the  segment, 
diaphragm  walls  must  be  placed  continuously,  centered  on  the  section  lines. 
Figure  2.1  shows  a  plan-view  cross  section  of  the  Braddock  Dam.  Steel  stiffeners 
may  be  placed  between  large  wall  panels  to  provide  support  for  external  pressure 
loads. 


2.2  Diaphragm  Panel  Construction 

The  diaphragm  walls  are  precast  flat  with  reinforcing  steel  protruding  to 
develop  a  connection  at  the  joints.  They  are  arranged  and  supported  in  their 
designed  location  before  the  bottom  slab  and  the  joints  are  cast  in  place.  Some 


2 


Chapter  2  In-the-Wet  Construction  Approach 


Figure  2.1 .  Plan  view  of  Braddock  Dam 


taper  is  required  at  most  joints  to  provide  a  stiffer,  stronger  section  where  the 
moment  and  shear  loads  are  greatest  and  to  provide  adequate  space  for  the 
confining  reinforcement.  This  design  should  be  standardized  as  much  as  possible 
in  order  to  reuse  formwork.  Figure  2.2  shows  a  typical  vertical  joint  placement. 


2.3  Critical  Load  Cases 

Critical  loads  for  most  components  of  the  structure  occur  when  the  fluid 
levels  in  adjacent  cells  are  substantially  different,  resulting  in  a  large  differential 
hydrostatic  pressure  on  the  concrete  slab.  The  load  cases  must  be  developed  in 
concert  with  the  construction  planner  in  order  to  identify  all  significant  planned 
and  accidental  events.  Each  concrete  slab  must  be  designed  for  both  positive  and 
negative  moments  along  both  slab  axes. 

Critical  loads  cases  defined  for  the  Braddock  Dam  were  these: 

a.  Float-out. 

b.  Outfitting  yard. 
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c.  Accidental  flood  event. 

d.  Touch-down. 

e.  Flood  event. 
f  Dead  load. 

g.  Pier  wall  in-fill. 

h.  Dewater  bay. 

i.  Second-stage  casting. 


2.4  Boundary  Conditions 

Four  boundary  conditions  must  be  modeled  in  the  typical  float-in  structure: 
continuous  elastic  foundation,  floating,  set-down  shafts,  and  drilled  piers.  One  or, 
in  some  cases,  two  or  three  of  these  conditions  are  in  effect  at  any  given  time. 
Set-down  shafts  and  drilled  piers  are  similar  in  their  design  and  function  with  the 
exception  that  only  set-down  shafts  are  engaged  when  the  structure  is  initially 
lowered  to  the  bottom  of  the  river. 

The  continuous  elastic  foundation  may  exist  during  initial  fabrication  in  the 
dry  dock  or  after  the  segment  is  placed  on  the  set-down  shafts  in  the  river  and 
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grouted.  Standard  finite  element  approaches  are  available  for  shell  elements  on 
an  elastic  foundation  (Chandrupatla  and  Belegundu  1997).  The  approximate 
stiffness  of  the  soil  or  grout  in  place  is  a  required  modeling  parameter.  The  force 
applied  to  a  structure  resting  on  an  elastic  foundation  is  proportional  to  the 
deflection  of  the  bottom  of  the  structure;  thus,  differential  settlement  is  expected. 
Typically,  the  upstream  end  of  a  navigation  structure  is  deeper,  so  it  will  have  a 
greater  weight  per  square  foot  of  plan  area. 

The  segment  will  be  floating  from  the  time  the  dry  dock  is  flooded  until  it  is 
placed  on  the  set-down  shafts.  The  pressure  applied  to  the  outside  of  the  structure 
is  a  function  of  the  draft  (the  depth  of  the  bottom  of  the  structure  when  floating). 
The  draft  is  based  on  the  buoyancy  of  the  structure  and  can  be  calculated  as 

D  =  Total  weiSht  (2 1) 

Plan  area  x  Water  density 

Pressure  is  exerted  on  all  submerged  faces  of  the  structure.  The  magnitude  of 
the  pressure  is  proportional  to  the  depth.  A  substantial  horizontal  pressure  can  be 
expected  on  the  outside  walls  of  the  structure,  which  is  not  present  when  the 
structure  is  supported  on  an  elastic  foundation.  A  finite  element  model  of  a 
floating  structure  should  be  in  equilibrium  if  the  pressure  is  input  exactly; 
however,  minimal  restraint  is  required  and  is  applied  by  including  springs  at  the 
four  comers  of  the  model.  The  construction  planner  must  ballast  the  structure  to 
account  for  higher  weight  per  square  foot  on  the  upstream  end  of  the  structure  so 
that  the  segment  floats  in  a  level  position. 

The  structure  is  initially  set  on  a  limited  number  of  drilled  piers.  These  are 
identified  as  set-down  drilled  shafts.  Additional  drilled  piers  are  distributed  under 
the  segment.  Precast  holes  in  the  bottom  slab  are  placed  over  drilled  piers.  These 
provide  lateral  and  vertical  support  after  the  hole  is  grouted.  These  shafts  in  a 
typical  structure  are  reinforced  concrete  on  the  order  of  6  ft1  in  diameter.  The 
shafts  can  be  modeled  as  springs  with  one  stiffness  value  in  the  vertical  direction 
and  another  in  the  horizontal  direction.  The  horizontal  stiffness  of  a  round  pier  is 
the  usually  the  same  in  any  direction.  The  vertical  stiffness  is  a  function  of  the 
stiffness  of  the  pier  itself,  the  foundation  material  at  the  bottom  of  the  shaft,  and 
any  support  gained  from  skin  friction.  The  translational  stiffness  is  controlled  by 
the  flexural  stiffness  of  the  reinforced  concrete  shaft  and  support  from  the 
surrounding  soil.  Methods  for  calculating  these  values  are  outlined  by  the  CASE 
Task  Group  on  Pile  Foundations  (Office,  Chief  of  Engineers  (OCE)  1983). 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  v. 
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3  Modeling  Considerations 


3.1  Assumptions 

The  structural  designer  must  know  the  overall  dimensions  of  the  segment,  the 
location  and  width  of  piers,  and  the  shape  of  the  individual  spillways  between 
piers.  The  designer  may  have  a  preliminary  layout  for  the  internal  diaphragm 
walls.  These  walls  are  required  to  support  the  external  slab  in  order  to  produce  an 
efficient  structure.  Spacing  between  diaphragm  walls  is  usually  fairly  uniform 
across  a  segment.  Diaphragms  should  be  spaced  to  break  the  slabs  into  sections 
that  are  nearly  square,  in  order  to  produce  two-way  slab  action.  Individual  walls 
may  be  placed  to  coincide  with  changes  in  the  geometry  of  spillways.  Diaphragm 
walls  are  closely  spaced  where  the  segment  rests  on  drilled  shafts  to  distribute  the 
concentrated  loads.  The  side  walls  of  a  pier  will  be  aligned  with  diaphragm  walls. 


3.2  Spillway  Shapes 


Spillway  layout  requires  the  selection  of  a  template  of  a  standard  spillway 
shape.  The  standard  spillway  shapes  are  based  on  typical  Corps  of  Engineers 
designs.  The  template  approach  assumes  that,  for  example,  the  outline  of  the 
typical  fixed-crest  weir  is  composed  of  the  same  series  of  geometrical  entities- 
flat  bottom  with  flat  ends  having  a  radius  at  the  comers  and  at  transitions 
between  flat  sections  on  the  top  surface.  The  template  fits  the  geometry  of  the 
segment  to  the  overall  dimensions  of  the  structure  by  extending  the  length  of  the 
level  area  on  the  downstream  end  of  the  spillway  to  the  predetermined  length  of 
the  segment.  A  minimum  number  of  changeable  parameters  are  used  to  describe 
the  geometry.  The  other  dimensions  of  the  spillway  are  calculated  using  basic 
geometrical  relationships  without  changing  the  essential  shape  of  the  spillway. 
Figure  3.1  shows  the  parameters  required  to  describe  the  shape  of  a  typical  fixed- 
crest  weir. 

The  spillway  shape  can  change  only  at  a  pier.  The  standard  design  approach 
for  a  pier  placed  where  the  spillway  shape  changes  is  to  extend  both  spillways 
into  the  pier.  For  example,  the  Braddock  Dam  segment  shown  in  Figure  2.1  has 
three  spillway  shapes.  The  pier  between  section  lines  E  and  F  transitions  from  the 
fixed-crest  weir  to  the  Water  quality  gate  bay.  This  can  create  complex  geometry 
for  modeling,  especially  when  the  top  surface  of  the  two  shapes  intersects  at  one 
or  more  points. 
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Figure  3.1 .  Fixed-crest  weir  template 


3.3  Tapers 

Tapers  are  typically  specified  at  points  where  slabs  intersect  at  angles  near 
90  deg.  They  are  used  to  stiffen  the  slab  where  the  highest  moment  and  shear 
loads  occur.  The  superelement  formulation  requires  at  least  a  minimal  taper  at  all 
comers  in  order  to  map  superelement  nodes  to  unique  shell  nodes.  Tapers  are 
defined  by  the  terms  “rise”  and  “run.”  These  terms  are  loosely  applied  because 
they  actually  define  dimensions  relative  to  surfaces  and  center  lines,  as  shown  in 
Figure  3.2.  The  location  of  the  taper  focal  point  (TFP)  is  defined  by  these  rise 
and  run  parameters.  The  run  is  the  distance  from  the  center  of  the  diaphragm  wall 
to  the  TFP,  and  the  rise  is  the  distance  normal  to  the  outside  surface  of  the 
spillway.  The  toe  of  the  taper  is  on  the  line  normal  to  the  outside  surface  through 
the  TFP,  as  shown.  The  rise  is  measured  parallel  to  that  line. 

To  enforce  uniformity  in  the  segment  layout  and  finite  element  model,  the 
rise  and  run  values  are  constrained.  Adjacent  tapers  must  be  consistent;  that  is, 
the  run  must  be  the  same  for  both  the  top  and  bottom  taper  on  one  side  of  a 
diaphragm  wall.  The  run  can  be  different  on  the  left  and  right  side  of  a  diaphragm 
wall.  Only  one  rise  value  is  permitted  along  the  bottom.  A  different  value  may  be 
applied  to  tapers  on  the  top. 
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3.4  Piers 

Piers  extend  above  the  spillway  surface  to  some  constant  elevation.  The 
walls  coincide  with  adjacent  diaphragm  walls  along  lettered  (N-S)  section  lines 
(Figure  2.1).  The  ends  are  enclosed.  The  upstream  end  is  usually  curved  for 
smoother  flow.  The  downstream  end  may  be  stepped  to  hold  temporary  steel 
gates  used  to  float  the  segment.  The  ends  of  pier  walls  are  modeled  as  flat  walls 
to  simplify  the  model. 

When  the  designer  designates  the  location  of  a  pier  wall  along  a  section  line, 
the  finite  element  model  is  divided  at  that  line.  The  shells  and  superelements  that 
would  normally  be  generated  along  the  section  line  are  divided  into  two 
appropriate  elements.  Figure  3.3  shows  the  finite  element  configuration  for  a 
simple  model  with  a  pier  wall.  Figure  3.3a  is  an  isometric  view  of  the  segment 
layout.  The  view  from  the  downstream  end  (Figure  3.3b)  shows  the  shell 
elements  in  the  model.  The  finite  element  model  is  divided  into  three  parts  (A,  B, 
and  C).  Figure  3.3c  is  the  end  view  of  the  model  with  the  superelements.  The 
boundaries  are  not  as  clear  in  this  sketch,  but  each  superelement  is  actually  two 
separate  elements. 

The  resulting  finite  element  model  for  the  segment  begins  as  several  separate 
segments  that  must  then  be  tied  together  in  the  finite  element  model.  The  tying 
operation  involves  identifying  all  shell  nodes  next  to  the  interface  and 
designating  one  set  as  masters  and  the  other  as  slaves.  Master  nodes  and  elements 
are  those  on  the  local  west  side  (-Z)  of  the  center  line.  Slaves  are  on  the  local  east 
side  (+Z).  Since  the  spillway  geometry  may  be  different  on  the  two  sides  at  the 
pier  wall,  these  nodes  will  generally  not  line  up.  However,  they  may  fall  on 
coincident  vertical  lines  since  the  diaphragm  locations  will  be  consistent.  The 
slave  node  degrees  of  freedom  are  mapped  to  the  master  elements  by 
(a)  determining  which  element  would  enclose  a  node  if  it  were  translated  to  the 
plane  of  the  master  elements,  (b)  calculating  the  local  coordinates  of  the  node 
location,  and  (c)  developing  a  transformation  matrix  to  map  slave  node 
displacements  to  the  master  element. 
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a.  Isometric  view  of  segment  layout 

Figure  3.3.  Finite  element  configuration  for  a  simple  model  with  pier  wall 
(Sheet  1  of  3) 
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b.  View  from  south  end  showing  three  submodels 
Figure  3.3.  (Sheet  2  of  3) 
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Model  A  Model  B  Model  C 


c.  View  from  local  south  end  (X)  with  solid  joint  elements 
Figure  3.3.  (Sheet  3  of  3) 


Figure  3.4  shows  a  slave  node  located  within  the  master  element  with  node 
numbers  1-8.  The  Z  coordinate  of  the  slave  node  is  one  half  the  thickness  of  the 
concrete  slab  from  that  of  the  master  element.  The  local  coordinates  in  the  master 
element,  £  and  t|,  are  found  by  inverting  the  shape  functions  for  the  shell  element 
iteratively. 

The  transformation  to  calculate  the  deflection  at  the  slave  node,  given  the 
deflections  of  the  nodes  in  the  master  element,  is  developed  in  two  steps.  The 
first  is  to  find  the  deflection  of  a  point  on  the  midplane  of  the  master  element 
adjacent  to  the  slave  node  using  die  standard  shape  function  relationship  (see 
Section  4. 1).  The  deflection  of  the  slave  node  is  then  mapped  to  that  point  by 
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where 


(3.1) 


10  0  0 
t 

0  1  0  -i- 

2 

0  0  10 
0  0  0  1 

0  0  0  0 


where  subscript  5  refers  to  slave  degrees  of  freedom,  and  subscript  m  refers  to 
masters.  The  thickness  of  the  slab  is  t.  The  rotations  at  the  two  points  are  equal, 
to  enforce  the  assumption  that  plane  sections  remain  plane  in  a  shell  element. 


Figure  3.4.  Slave  node  is  mapped  in  to  master  element 


Slave  nodes  are  eliminated  from  the  finite  element  solution,  and  the  stiffness 
matrix  for  elements  containing  slave  nodes  is  transformed  to  align  with  the 
master  degrees  of  freedom  by  developing  a  transformation  matrix  for  all  degrees 
of  freedom  in  elements  with  slave  nodes  using  the  transformations  described 
above.  The  deflections  at  each  slave  nodes  are  a  function  of  the  deflections  at  up 
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to  eight  master  nodes.  The  transformation  matrix,  [Teiem],  is  assembled  from 
[Tnode]  submatrices  weighted  by  the  shape  functions  relative  to  the  master  node. 

kl-fcJkl  <3-2> 

The  stiffness  matrix  for  the  slave  element  is  then  calculated  from 


[/^master  ]  -  1 Teiem  F  slave  ]  I Teiem  ] 


(3.3) 


3.5  Loading 

Loads  on  the  structure  include  gravity  loads  imposed  as  body  forces  on  the 
elements,  external  pressures  applied  inside  the  cells  or  on  the  outside  walls  of  the 
structure,  and  concentrated  loads  that  can  be  applied  to  nodes.  Body  forces  are 
calculated  at  the  integration  points  and  distributed  to  the  nodes  by 

{P„}  =  2\){f}Tbr\j\d^  (3.4) 

-1-1 

Hydrostatic  pressure  loads  are  defined  for  various  load  cases  by  designating 
the  fluid  level  in  each  cell  or  external  face  relative  to  the  bottom  of  the  structure 
(Y  =  0).  The  fluid  may  be  water  or  tremie  concrete.  When  a  hydrostatic  pressure 
is  imposed  on  a  shell  element,  the  depth  of  the  integration  point  is  calculated 
relative  to  the  top  of  the  fluid,  and  the  pressure  is  calculated  by  multiplying  the 
depth  by  the  density  of  the  fluid.  This  pressure  is  then  distributed  to  the  nodes 
using  the  shape  functions  by 

{rj-A'fyY  faWWi  <3-5> 

-1-1 

The  total  pressure  is  distributed  to  the  x,  y,  and  z  directions  by  multiplying  by  the 
direction  cosines. 

This  approach  assumes  that  the  pressure  varies  linearly  in  the  vertical 
direction.  This,  of  course,  is  not  exact  when  the  top  of  the  fluid  is  in  the  middle 
of  an  element.  In  this  case,  the  pressures  at  the  integration  points  are  adjusted  to 
impose  the  same  total  force  but  distribute  it  across  the  entire  element.  The 
moment  will  tend  to  be  greater  in  this  case,  but  the  errors  decrease  as  the  finite 
element  mesh  is  refined. 
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4  Shell  Element  Formulation 


An  eight-node,  40-degrces  of  freedom  (DOF),  biquadratic  shell  element  is 
used  to  model  the  concrete  slabs.  The  formulation  follows  Weaver  and  Johnston 
(1984).  The  global  finite  element  model  uses  six  DOF  per  node.  Stiffness  at  the 
translation  degrees  of  freedom  is  calculated  directly  in  the  global  coordinate 
system  using  a  transformed  materials  property  matrix.  The  local  rotational  DOF, 
a  and  J5,  are  transformed  based  on  the  direction  cosines  at  each  node. 


4.1  Geometry 

Figure  4.1  shows  the  node  numbering  and  natural  (0  t|,  0  coordinate  system 
for  the  element.  The  two  rotational  DOF,  a  and  P,  are  aligned  with  the  local  x’ 
and  y’  directions  of  the  element  at  any  point.  The  z’  direction  is  normal  to  the 
midplane  of  the  element.  The  x’  and  y’  directions  are  perpendicular  to  z’  and  to 
each  other  at  any  point. 

The  local  coordinate  directions  are  stored  as  direction  cosines  at  the  nodes 
and  integration  points.  These  are  used  throughout  the  finite  element  formulation. 
The  direction  cosines  are  designated  as 


1 2  (  3 

m2  mi 
n2  «3 


(4.1) 


where  mi;  and  nj  are  the  x,  y,  and  z  components  of  the  unit  vector  Vj.  In  this 
formulation,  V3  is  normal  to  the  midplane  of  the  shell,  and  V,  and  V2  are  aligned 
with  the  midplane  normal  to  V3  and  to  each  other. 

The  location  of  any  point  in  the  element  can  be  calculated  from  the  natural 
coordinates  as 


X 

8 

V 

+i/4 

'  /*  ‘ 

y 

y, 

”3  i 

z 

i=l 

.zi. 

/=!  * 

3/_ 

(4.2) 


14 


Chapter  4  Shell  Element  Formulation 


Figure  4.1.  Shell  element 


where  the  shape  functions  f  are  functions  only  of  the  in-plane  coordinates. 
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Similarly,  displacements  are 
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4.2  Stiffness 


The  derivation  of  the  stiffness  matrix  follows  typical  finite  element 
procedures.  The  Jacobian  matrix  for  a  three-dimensional  element  is 

Z’t 

[J]=  x>n  y>n  z>n  I 

x,(-  y,^  zi( 

The  3x3  Jacobian  matrix  is  inverted  to  find  the  derivatives  of  the  natural 
coordinates  in  terms  of  the  global  coordinate 


U’x  *hx  C*x 

K*]=  Thy 

TJ’z  C’z 


The  strain-displacement  matrix  is  derived  by  defining  individual  strain 
components  in  terms  of  the  partial  derivatives  of  the  displacements. 


U,y  +  v,x 

v.,  +  ^ 
w.  +u. 


A  6x5  submatrix  of  the  6x40  strain-displacement  matrix  is  then 


o,  0  0 

0  b,  0 

0  0  c, 


~  eim2i 


r5  ]=  u  u  c-  -SKi  Sinu 

bi  ai  0  -e,*2  i-dim2i  ei£u+dimu 

0  c,  bt  -  g.m2i  -  etnv  g,mu+e,nu 

-Ci  0  °i  ~din2,  ~g^2i  d,nu+g/-U 


0  =  1,2,...,  8) 
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where 


=  J\\fi4  +  Jufifl 

<7,=|(<  +  4/l) 

=  *^21  fi,4  +  Jllfifl 

=  ^  3\f ij;  +  J  32  fij 

a  =|(^+4/,) 

The  stress-strain  relationship  (Equation  4.8)  is  a  function  of  the  Young’s 
modulus  (E)  and  shear  modulus  (G)  of  the  material.  The  third  row  and  column 
are  zero  in  the  local  coordinate  system  because  there  is  no  stress  through  the 
thickness  (plane  stress  assumption). 


E,,  0  0  0  0 

£,v  E,y.  0  0  0  0 

0  0  0  0  0  0 

0  0  0  Gxy  0  0 

G 

0  0  0  0  0 

1.2 

C* 

0  0  0  0  0  — ^ 

1.2  . 

This  relationship  is  transformed  to  the  global  coordinate  system  by 

[£]=[rj[£'][rj  (4.9) 


where 


l\  wf  lxn\  777,77,  77,/, 

l\  nt,  «z4 
l]  rr%  /J32  /3m,  Wh  fhh 

2//2  2m, 2/7,^  +l2n\  n^r^+mft  r\l2  +n2ll 

2/24  2^7773  2^773  +4m,  mfo+mfo  nj3  +77j4 

2/3/]  2777,777,  277,77,  Ifa  +47^  777,77,  +777,773  nJi+rhh 


(4.10) 
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Finally,  the  stiffness  matrix  is  calculated  by  integrating  [B]T[E][B]  over  the 
volume  of  the  element.  This  operation  must  be  performed  by  numerical 
integration  by  summing  the  values  at  the  integration  points.  Four  integration 
points  are  used  in  the  current  formulation.  These  are  on  the  midplane  (£  =  0).  The 
in-plane  natural  coordinates  of  the  four  points  are 


4,  =- 


V3 


£  =  — 
V3 


7j  =  V3 


(4.11) 


,  _  1  1 
*‘"V3  ,<  =  Ji 


The  strain-displacement  matrix  is  divided  into  two  components  to  separate 
out  the  terms  containing  the  £  term: 


The  stiffness  matrix  as  calculated  by  Equation  4.13: 

m=z[2 


K’->h 


(4.13) 


4.3  Load  Boundary  Conditions 


4.3.1  Body  forces 

Body  forces  are  due  to  gravity  in  the  negative  Y  direction.  Equivalent  nodal 
loads  due  to  gravity  body  forces  with  a  material  density,  by,  are 


=  *L 
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4(^,7,.) 

4(^,7,.) 


by\AWi 


(4.14) 


18 


Chapter  4  Shell  Element  Formulation 


4.3.2  Pressure 


When  pressure  is  applied  to  an  element  face,  the  equivalent  nodal  forces  are 
determined  based  on  the  pressure  function  at  the  integration  points.  The  total 
force  at  the  node  is  calculated  as 
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mnv>) 

J73  jra 
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-2 

/=! 

Pey 

/•«,.*,) 

Piy 

;,) 

kJ 

/•(£,. 7, X 

P(yMw, 


(4.15) 


Pressure  components  in  the  x,  y,  and  z  directions  are  determined  by 
multiplying  the  P(yj)  by  the  direction  cosines  at  the  integration  point. 

If  an  element  is  only  partially  submerged,  this  approach  would  accurately 
model  the  linearly  varying  positive  pressure,  but  it  would  effectively  apply  a 
negative  pressure  to  the  element  above  the  fluid.  A  simple  approximation  is 
employed  to  adjust  the  nodal  loads  to  duplicate  the  total  pressure  force  on  the 
element  without  creating  negative  pressures.  The  node  with  the  greatest  elevation 
is  determined  for  the  element.  The  fluid  level  is  assumed  to  be  at  this  point,  and 
the  density  of  the  fluid  is  reduced  to  exert  the  same  force. 

For  example,  a  rectangular  shell  element  is  20  ft  wide.  The  elevation  at  the 
bottom  of  the  element  is  10  ft,  and  the  elevation  is  25  ft  at  the  top.  The  fluid  level 
is  20  ft  of  water  with  a  density  of  62.4  pcf.  The  total  pressure  force  on  the 
element  is 


P  =  (20  ft)  (20  ft  -10  ft)  (10  ft  x  62.4  lb/ft3)/2  =  62,400  lb 


Multiplying  the  density  by 


20’-10'Y 

25’— 10' J 


4 

=  —  reduces  the  density  to  27.7  pcf  and 


imposes  the  same  force.  The  force  distribution  will  be  somewhat  inaccurate,  but 
the  error  will  decrease  as  the  mesh  is  refined. 
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5  Superelement  Formulations 


The  DBA  approach  uses  superelements  to  accurately  model  the  geometry  at 
joints  between  concrete  slabs.  When  multiple  shell  elements  meet  at  a  point,  the 
implied  geometry  does  not  duplicate  the  actual  geometry  of  a  tapered  joint.  The 
superelement  approach  allows  the  designer  to  study  the  effect  of  various  taper 
geometries  and  to  evaluate  three-dimensional  (3-D)  states  of  stress  in  the  joints. 
Standard  element  shapes  are  defined  parametrically,  so  actual  structural 
dimensions  are  used  to  generate  the  refined  mesh  and  develop  the  superelement 
stiffness  matrices.  The  superelement  formulation  involves  the  creation  of  a 
refined  finite  element  model  consisting  of  solid  hexahedral  and  prism  elements. 
Degrees  of  freedom  at  nodes  that  do  not  interface  with  other  elements  in  the 
structure  are  removed  using  static  condensation.  These  are  designated  as  slave 
nodes.  The  remaining  master  nodes  are  then  constrained  to  mate  with  the  6-DOF 
shell  nodes  in  the  rest  of  the  model.  These  procedures  are  computationally 
intensive,  so  a  novel  approach  has  been  developed  to  synthesize  the  superelement 
stiffness  matrix  and  the  transformation  matrix  that  maps  the  response  of  slave 
nodes  to  the  master  nodes.  These  matrices  for  an  arbitrary  element  configuration 
are  assumed  to  be  linear  combinations  of  a  set  of  baseline  matrices.  The 
coefficients  used  in  the  linear  combination  are  some  function  of  the  physical 
parameters  that  define  the  superelement,  for  example,  length  and  modulus.  These 
functions  are  determined  by  training  neural  networks  using  actual  stiffness 
matrices. 


5.1  Standard  Elements 

Six  standard  superelements  have  been  developed.  These  are  categorized  as 
point  (PI,  P2,  and  P4)  and  line  (L2,  L3,  and  L4)  elements.  Figure  5.1  shows  a 
simple  structure  with  examples  of  each  element  type  labeled. 

Point  elements  join  the  comers  of  shells.  Three  types  are  available  (PI,  P2, 
and  P4)  with  three,  five,  and  eight  shell  nodes,  respectively,  depending  on  their 
location  in  the  model.  They  are  characterized  by  five  parameters:  depth  (D), 
width  (W),  height  (H),  rise  (S),  and  run  (N).  Figure  5.2  gives  details  of  the 
characteristics  of  each  element.  The  shell  nodes  are  numbered;  the  parameters  are 
defined  in  Figure  5.2c.  These  five  parameters  were  selected  to  adequately  model 
typical  joint  geometry.  Minor  approximations  are  made  in  some  cases;  for 
example,  the  run  (N)  may  vary  on  different  sides  of  an  element.  An  average  value 
is  used  to  limit  the  number  of  parameters  describing  the  superelement.  This 


20 


Chapter  5  Superelement  Formulations 


method  has  the  inherent  flexibility  to  add  parameters  in  order  to  provide  more 
modeling  options. 


Figure  5.1 .  Typical  superelement  locations 


Figure  5.3  shows  the  three  types  of  line  elements.  These  are  designated  L2, 
L3,  and  L4  with  six,  nine,  and  twelve  nodes,  respectively.  They  interface  with 
two,  three,  or  four  shell  elements  along  the  sides  and  point  elements  on  the  ends. 
The  midside  nodes  are  provided  for  compatibility  with  the  eight-node  shell 
elements.  The  line  elements  are  described  by  four  parameters:  length  (L), 
width  (W),  depth  (D),  and  rise  (S)  (as  labeled  in  Figure  5.3c).  Again,  the  number 
of  parameters  was  limited  to  facilitate  the  development  of  the  neural  networks  in 
the  initial  development  of  this  method.  In  this  case,  rise  and  run  are  constrained 
to  be  equal  on  all  sides  of  the  element. 


Chapter  5  Superelement  Formulations 


21 


\ 

s 

si 

s 

s 

s. 

S 

s. 

s 

s. 

s 

x 

S 

SB 

V 

s 

s 

S 

II 

\ 

§ 

s 

tt— 

s 

- 

s 

S| 

\. 

\ 

A 

s 

x: 

s 

■Bl 

_ Ji 

\ 

B 

a.  Point  superelement  PI 


b.  Point  superelement  P2 

Figure  5.2.  Point  superelements  (Continued) 
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c.  Point  superelement  P4 
Figure  5.2.  (Concluded) 


a.  Line  superelement  L2 

Figure  5.3.  Line  superelements  (Continued) 
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c.  Line  superelement  L4 
Figure  5.3.  (Concluded) 
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5.2  Stiffness  Matrix  Calculation 

The  calculation  of  the  superelement  matrices  begins  with  the  assembly  of  a 
conventional  stiffness  matrix  for  the  solid  model.  A  simple  preprocessor  is  used 
to  generate  the  finite  element  mesh  for  each  element  type.  A  typical  coarse  mesh 
is  plotted  in  Figures  5.2  and  5.3  above.  Standard  eight-node,  24-DOF  hexahedral 
elements  are  used.  Prism  elements  with  four  quadrilateral  and  two  triangular 
faces  are  created  as  a  degenerate  hexahedral  element  with  two  pairs  of  coincident 
nodes. 

After  assembling  the  stiffness  matrix  for  the  solid  model,  individual  nodes 
are  designated  as  slaves  and  masters.  Master  nodes  are  those  on  a  surface  that 
interfaces  with  a  shell  element.  All  other  nodes  are  slave  nodes.  A  Gaussian 
elimination  procedure  is  applied  to  the  stiffness  matrix  to  create  both  the 
condensed  stiffness  matrix  and  a  displacement  transformation  matrix  (T).  The 
procedure  is  documented  in  many  finite  element  references  (e.g.,  Weaver  and 
Johnston  1984,  Bathe  1982).  The  typical  mathematical  formulation  involves 
rearranging  the  degrees  of  freedom  and  partitioning  the  stiffness  matrix  into  slave 
and  master  submatrices. 


KSs 

K  MS 


Ksm 
K-MM  _ 


(5.1) 


The  [Tsm]  matrix  is  used  to  calculate  the  deflections  at  slave  nodes  given  the 
results  at  master  nodes. 

{us  }  =  I Tsm  K  }  ($-2) 


The  [TSm]  matrix  can  be  derived  from  the  submatrices  as 
\?sm  ] =  —  [^ss  1  [^\w] 


(5.3) 


The  condensed  stiffness  matrix  is  then 


(5.4) 


where 

~  [^SS  r  \Ksm  ]  /r  c 

.  M  J  ( 

and  [I]  is  the  identity  matrix. 

The  Gaussian  elimination  procedure  developed  to  create  superelements 
generates  these  matrices  without  performing  the  complete  inversion  of  [KSs],  but 
it  is  still  computationally  intensive. 
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Following  the  condensation,  the  stiffness  matrix  must  be  constrained  to 
comply  with  the  six-DOF  shell  nodes  in  the  global  finite  element  model. 

Figure  5.4  shows  the  master  nodes  on  the  typical  line-shell  element  interface.  The 
shell  nodes  always  coincide  with  a  master  node  in  the  line  element.  This  node, 
along  with  the  adjacent  nodes  on  the  top  and  bottom  surface  of  the  shell,  must 
remain  on  a  line  to  enforce  a  fundamental  assumption  in  the  development  of  the 
shell  element. 
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1  0  0 
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(5.6) 
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Other  nodes  along  the  interface  are  mapped  to  the  shell  nodes  using  a 
quadratic  shape  function.  A  transformation  matrix,  [Tc],  is  developed  to  relate  the 
displacement  of  the  master  nodes  to  the  shell  degrees  of  freedom.  The  final 
superelement  stiffness  matrix  with  shell  nodes  is  then  calculated  as 


[**.]=  (5.7) 

The  transformations  used  to  develop  the  superelement  are  also  used  for  result 
post-processing.  Displacements  for  the  shell  nodes  in  a  superelement  are 
extracted  from  the  global  finite  element  results.  A  transformation  matrix  is 
generated  to  convert  these  to  X,  Y,  and  Z  displacements  at  each  node  in  the 
original  solid  model.  This  matrix  is  formed  by  multiplying  the  [TSe]  and  [Tc] 
matrices  described  above. 

k^}=[r»][r.]k*}  (5-8) 

Strains  and  stresses  are  then  calculated  using  the  strain-displacement  [B] 
matrix  for  the  solid  elements  and  a  standard  3-D  material  property  matrix. 


5.3  Neural  Network  Synthesis 

The  development  of  superelement  matrices  is  computationally  intensive.  A 
typical  navigation  structure  model  will  have  on  the  order  of  1,000  superelements, 
so  faster  methods  are  desired.  However,  the  ability  to  accurately  model  the  true 
configuration  and  response  of  the  structure  must  be  maintained.  Methods  that 
allow  computationally  intensive  operations  to  be  performed  one  time  during 
program  development  and  make  the  results  available  in  a  design  tool  can  provide 
an  efficient  approach. 

The  typical  innovative  navigation  structure  can  be  accurately  modeled  using 
the  shell,  point,  and  line  elements  described  above.  The  range  of  sizes  and  shapes 
of  these  elements  is  somewhat  restricted  by  the  nature  of  the  design.  A  catalog  of 
elements  could  be  developed,  and  the  program  Could  select  the  most  appropriate 
element  at  run  time.  This  approach  would  provide  good  results  if  enough 
elements  are  available  in  the  catalog.  However,  it  would  be  tedious  and 
inefficient  to  manage  the  large  database  required,  and  the  accuracy  would  always 
depend  on  how  close  the  actual  element  is  to  a  cataloged  element. 

It  would  be  desirable  to  be  able  to  work  with  a  small  catalog  of  elements  and 
interpolate  between  them  to  form  a  superelement  that  accurately  models  the 
structure  geometry.  The  approach  developed  in  this  research  involved  creating  a 
set  of  N  “baseline”  matrices,  [KJ,  representing  the  geometry  of  several  elements 
throughout  the  parameter  space.  The  matrices  for  a  new  superelement  would  be 
created  as  a  linear  combination  of  the  stiffness  and  transformation  matrices  of  the 
baseline  matrices. 
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Even  with  dozens  of  baseline  matrices,  this  operation  is  much  faster  than 
conventional  methods. 

Appropriate  coefficients  must  be  calculated  based  on  the  parameters  for  the 
desired  superelement.  In  the  simple  case  of  a  one-parameter  model,  these  could 
be  found  by  calculating  the  value  of  each  coefficient  as  a  function  of  that 
parameter  for  several  elements  and  then  fitting  a  curve  through  those  points  using 
a  regression  analysis.  Figure  5.5  shows  a  plot  of  a  simple  application  of  this 
method  using  linear  regression.  Since  multiple  parameters  define  the 
superelements,  and  the  form  of  the  function  that  would  best  fit  the  data  is  not 
known,  a  neural  network  approach  was  selected  to  develop  this  relationship. 
Figure  5.6  is  a  schematic  diagram  of  a  back-propagation  neural  network  used  to 
develop  a  function  to  calculate  coefficients  from  the  geometric  parameters.  These 
methods  are  well  documented  (see,  for  example,  Fausett  1994;  Hagan,  Demuth, 
and  Beale  1996;  Masters  1993). 


Coefficent  Number  1 


Figure  5.5.  Linear  regression  to  estimate  a  coefficient  given  the  length  (sample) 


A  neural  network  is  usually  defined  as  an  electronic  circuit  with  many  simple 
processors  that  attempts  to  mimic  the  function  of  the  brain.  It  can  also  refer  to  a 
computer  algorithm  that  duplicates  the  function  of  that  circuit.  The  purpose  of  the 
network  is  to  gain  “knowledge”  through  experience.  That  knowledge  is  stored  as 
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Figure  5.6.  Schematic  diagram  of  a  neural  network 


synaptic  weights.  In  this  application,  the  neural  network  is  expected  to  leam  the 
values  of  the  coefficients  used  in  the  linear  combination  given  the  input 
parameters.  It  gains  experience  by  being  fed  many  sets  of  input  and  output  data. 

The  development  of  the  neural  network  model  is  a  two-step  process  after  the 
baseline  matrices  have  been  determined.  First,  the  “best”  coefficient  values  must 
be  calculated  for  several  superelements,  after  which  a  neural  network  must  be 
trained  using  these  data.  The  coefficient  values  are  calculated  using  a  least- 
squared  error  approach.  This  is  the  same  approach  used  in  linear  regression.  In 
linear  regression,  values  are  calculated  for  the  slope  and  y-intercept  of  the  line 
that  minimize  the  sum  of  the  squared  error  between  the  actual  data  and  the  line. 
In  this  development,  the  error  is  the  sum  of  the  relative  errors  between  each  term 
in  the  actual  matrix  and  the  corresponding  term  in  the  linear  combination, 


where 


M  =  number  of  DOF  in  the  superelement 
N  =  number  of  baseline  matrices 
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Kjk  =  terms  in  the  actual  stiffness  matrix 
Kijk  =  elements  in  jth  baseline  matrix. 

Partial  derivatives  of  the  error  function  are  taken  with  respect  to  each 
unknown  coefficient,  ci5  and  set  equal  to  zero.  This  set  of  simultaneous  equations 
is  then  solved  to  determine  the  coefficients.  The  accuracy  of  the  matrix  calculated 
by  the  linear  combination  is  improved  with  more  baseline  matrices.  The  target 
error  was  an  average  less  than  1  percent  with  a  maximum  error  less  than 
10  percent.  This  was  achieved  with  about  50  baseline  matrices. 

Each  line  in  the  network  is  assigned  a  weighting  function  that  is  applied  to 
the  value  of  the  input  (left)  node  in  order  to  calculate  the  excitation  at  the  next 
(right)  node.  One  or  more  layers  of  hidden  nodes  are  used  to  improve  the 
accuracy  of  the  model.  The  network  must  be  trained  in  order  to  determine  the 
correct  weighting  functions.  This  process  requires  training  data  derived  from 
actual  superelements.  The  training  data  are  input  parameter  values  paired  with 
the  coefficients  required  to  calculate  the  superelement  matrix  from  the  linear 
combination  of  baseline  matrices.  Several  hundred  sets  of  data  are  generated  to 
train  the  neural  network.  Initial  weighting  functions  are  arbitrarily  assigned  to  the 
network. 

The  first  set  of  input  parameters  is  run  through  the  network  to  calculate  an 
approximate  set  of  outputs.  The  error  between  the  approximate  outputs  and  the 
correct  outputs  is  calculated,  and  the  weighting  functions  are  adjusted  to  reduce 
this  error  in  the  back-propagation  phase.  The  sets  of  training  data  are  processed 
through  the  network  many  times  until  the  error  is  reduced  to  an  acceptable  level. 
The  weighting  functions  are  then  stored,  and  this  trained  neural  network  becomes 
the  function  used  to  calculate  the  coefficients  used  in  the  linear  combination  to 
determine  the  superelement  matrices. 
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6  Design  Considerations 


The  design  of  innovative  navigation  structures  is  based  on  the  requirements 
of  ACI  318-02.  Slabs  are  designed  using  the  procedures  for  two-way  slabs.  Body 
forces  are  considered  as  dead  load,  and  all  hydrostatic  pressures  and  concentrated 
loads  are  assumed  to  be  live  loads. 


6.1  Two-Way  Slab  Design 

The  ACI  318-02  provisions  for  flexure  design  of  two-way  slabs  as  they  are 
applied  for  DBA  are  summarized  in  Appendix  A.  When  determining  the 
minimum  slab  thickness  (see  Appendix  A,  Section  A. 6),  it  is  assumed  that  walls 
support  all  slabs.  In  this  case,  the  stiffness  of  the  supporting  “beam”  is  large,  so 
ctn,  is  assumed  to  be  greater  than  2  in  all  cases. 

Since  the  maximum  moment  is  expected  be  greater  in  the  short  direction,  the 
reinforcement  in  this  direction  of  the  slab  is  placed  at  the  outside  on  both  the  top 
and  bottom  with  the  required  clear  cover.  Reinforcement  in  the  long  direction  is 
on  the  inside,  resulting  in  a  smaller  depth  of  reinforcement.  The  reinforcement  in 
each  slab  is  based  on  the  worst-case  combination  of  shear,  moment,  and  thrust  at 
any  point  in  that  slab  for  any  specified  load  combination.  As  a  minimum, 
reinforcing  bars  are  placed  at  12-in.  spacing  on  both  faces  in  all  slabs  to  provide 
the  minimum  reinforcement  ratio  (0.0018  for  Grade  60).  Slabs  up  to  9  in.  thick 
require  No.  3  bars.  Slabs  between  9  and  16  in.  use  No.  4  bars  at  12  in.  The  size 
and  spacing  of  bars  is  modified  based  on  structural  requirements. 

The  designer  provides  the  initial  thickness  of  the  slab.  Given  the  moment  in 
the  short  direction  of  the  slab,  designs  are  developed  for  each  available  bar  size, 
since  the  depth  of  reinforcing  is  a  function  of  bar  size.  Assuming  that  the  largest 
allowable  bar  size  should  be  used  to  minimize  the  number  of  bars  in  construction, 
the  design  with  the  largest  bar  size  that  satisfies  the  requirements  for  maximum 
bar  spacing  is  selected.  This  bar  diameter  is  then  added  to  the  clear  cover  for  the 
design  of  reinforcing  in  the  long  direction  of  the  slab,  and  the  procedure  is 
repeated  for  moments  in  that  direction. 

Stresses  caused  by  thrust  loads  will  tend  to  be  much  lower  than  those  caused 
by  flexure.  ACI  requirements  for  minimum  reinforcement  in  walls  are  checked 
for  all  slabs,  but  it  is  unlikely  that  these  will  govern.  A  simple  procedure  was 
developed  to  superimpose  the  two  cases.  When  the  thrust  is  in  compression,  the 
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compressive  strength  of  the  concrete  used  in  the  flexural  design  calculation  is 
reduced  by  the  compressive  stress,  oc,  factored  based  on  wall  design 
requirements. 


'  c,  reduced 


=/;- 


(0.7X0.75) 


When  the  thrust  load  is  in  tension,  the  reinforcement  is  increased  in  both  the 
top  and  bottom  of  the  slab  to  resist  the  additional  force. 

Checks  are  then  performed  to  ensure  that  the  reinforcement  ratio  is  less  than 
75  percent  of  the  balanced  steel  ratio  in  flexure.  If  the  steel  ratio  is  too  high,  the 
thickness  of  the  slab  must  be  increased.  Since  the  design  objective  is  a 
lightweight  structure,  the  reinforcement  ratio  should  be  near  the  maximum. 

The  shear  capacity  of  the  slab  is  a  function  of  thickness.  It  is  assumed  that 
reinforcing  will  be  used  only  for  shear  in  local  areas.  This  shear  capacity  is 
factored  and  compared  with  the  shear  at  every  analysis  point. 

K  =  (0.85)2  tfbd  (6.2; 

The  factor  of  safety  (F.S.)  must  be  greater  than  1  at  all  points  in  the  slab. 


F.S.  = 


where  Vfe  is  the  shear  load  per  foot  determined  from  the  finite  element  analysis, 
and  the  width,  b,  is  12  in.  This  procedure  takes  advantage  of  the  finite  element 
results  to  develop  a  simpler  algorithm,  but  it  is  consistent  with  the  assumptions  in 
ACI318. 


6.2  Superelement  Analysis 

The  use  of  superelements  in  the  joints  of  the  structure  permits  the  designer  to 
study  the  stresses  in  more  detail.  In  the  typical  design,  the  required  reinforcing 
for  the  slabs  is  extended  into  the  joints,  and  some  additional  reinforcing  is  added 
to  the  tapered  regions.  DBA  methods  have  not  been  developed  to  detail  joint 
reinforcing. 


6.3  Controlling  Case  Determination 

The  design  of  a  slab  is  based  on  the  reinforcing  requirements  at  the  critical 
point  in  the  slab  for  all  load  cases.  Each  slab  must  be  checked  for  all  load  cases. 
The  case  requiring  the  most  reinforcement  in  each  of  the  four  layers  will  govern 
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6.4  Automatic  Design  Modification 

An  initial  slab  thickness  must  be  proposed  by  the  designer  in  order  to  analyze 
the  structure.  The  initial  design  is  based  on  these  analysis  results.  These  results 
will  suggest  changes  to  the  design  in  some  situations.  One  goal  in  the  design  is  to 
produce  the  lightest  structure  that  meets  all  requirements.  A  second  goal  is  to 
produce  a  constructible  design  that  allows  for  reuse  of  formwork  and  reinforcing 
schedules  that  are  easy  to  fabricate.  The  automatic  redesign  procedures 
developed  for  DBA  are  based  on  these  two  goals. 

Designs  are  screened  to  identify  slabs  that  may  be  overdesigned.  These  are 
slabs  with  a  high  factor  of  safety  for  shear  and  a  low  reinforcement  ratio.  The 
thickness  of  these  slabs  is  reduced  in  1-in.  increments  and  redesigned  based  on 
the  shear,  moment,  and  thrust  values  from  the  initial  analysis.  Although  this  is 
not  exact,  small  changes  in  thickness  will  have  little  effect  on  these  results.  If  the 
design  for  this  thinner  section  is  acceptable,  the  slab  thickness  is  changed  in  the 
design.  The  thickness  is  further  reduced  and  the  redesign  calculations  are 
repeated  until  an  unacceptable  design  is  produced. 

Similarly,  automatic  redesign  procedures  increase  the  thickness  of  slabs  that 
are  underdesigned.  The  shear  factor  of  safety  is  less  than  1  in  these  slabs,  the 
steel  reinforcement  ratio  is  greater  than  the  allowable,  the  thickness  is  insufficient 
to  construct  the  slab  given  the  required  reinforcement,  or  the  thickness  is  less 
than  the  minimum  allowed  for  serviceability  by  ACI 318. 
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7  Shear,  Moment,  and  Thrust 
Calculation 


Design  and  analysis  methods  for  concrete  structures  are  based  on  the  shear, 
moment,  and  thrust  (axial  forces)  actions  on  the  beam,  wall,  or  slab  section. 
Various  methods  are  available  to  approximate  these  actions.  The  finite  element 
method  is  used  in  this  development. 

The  design  of  innovative  navigation  structures  is  concerned  with  ensuring 
that  the  thickness  and  reinforcement  in  the  slab  sections  is  adequate  to  resist  the 
applied  loads.  These  sections  are  treated  as  two-way  slabs  in  the  design.  Stresses 
can  be  calculated  in  the  joints,  but  design  criteria  are  not  as  well  established 
when  the  overall  response  is  more  complex  than  simple  bending  actions.  The 
superelement  approach  developed  in  this  research  allows  the  designer  to  evaluate 
these  stresses.  Capabilities  are  also  included  to  evaluate  overall  bending  action  in 
a  cross  section  of  the  structure. 


7.1  Shell  Elements 

Methods  for  calculating  shear,  moment,  and  thrust  in  shell  elements  are  well 
established.  Because  the  individual  stress  components  will  be  used  in  later  cross- 
section  plots,  these  were  also  calculated.  Actions  in  each  shell  element  are 
evaluated  by  extracting  the  nodal  displacements  from  the  global  finite  element 
results.  The  global  rotations  (etc,  Pg,  and  yo)  are  transformed  back  to  the  local  (a 
and  p  only)  system.  The  strain-displacement  [B]  matrix  (Equation  4.12)  is  then 
evaluated  at  each  integration  point  in  the  local  coordinate  system  and  used  to 
calculate  the  six  local  strain  components  at  the  bottom  (£  =  -1)  and  top  (£  =  +1) 
of  the  element.  Stresses  are  then  determined  from  the  strains  and  used  to 
calculate  the  moment,  shear,  and  thrust  at  the  integration  points.  Since  the  stress 
results  are  aligned  with  the  local  shell  element  coordinates,  this  calculation  is 
straightforward. 

Defining  the  element  coordinates  as  1  and  2,  the  shear,  moment,  and  thrust  in 
the  1  direction  can  be  determined  from  the  extensional  stress  in  that  direction  and 
the  shear  stress  in  the  13  plane.  Actions  in  the  2  direction  are  derived  from 
stresses  in  that  direction  (Figure  7.1).  Assuming  that  the  stress  varies  linearly 
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through  the  thickness,  the  thrust  component  is  a  function  of  the  average  of  the 
stresses  on  the  top  and  bottom.  This  average  stress  is  multiplied  by  the  thickness 
to  determine  the  thrust  per  unit  length.  If  the  thrust  component  is  removed  from 
the  stress  distribution,  die  resulting  stress  is  pure  bending.  The  moment  per  unit 
length  is  then  simply  a  function  of  the  maximum  stress.  Finally,  shear  per  unit 
length  is  computed  from  the  appropriate  shear  stress  component  by  multiplying 
by  the  thickness.  The  actual,  parabolic  shear  stress  distribution  cannot  be 
calculated  in  a  standard  shell  element  because  the  derivation  is  based  on  the 
assumption  that  plane  sections  remain  plane;  however,  the  calculation  of  total 
shear  is  accurate.  Figure  7.2  shows  the  derivation  of  thrust  and  moment. 

Once  shear,  moment,  and  thrust  are  calculated  at  the  integration  points,  they 
are  extrapolated  to  the  nodes  based  on  linear  shape  functions.  Essentially,  a  four- 
node,  bilinear  element  is  constructed  with  the  integration  points  as  the  comers. 
The  standard  shape  functions  for  this  element  are  used  to  extrapolate  to  the  nodes 
in  the  shell  element.  For  example,  Node  II  in  Figure  7.3  is  at 

r  =  V3  and  s  =  — \[3  .  The  actions  at  Node  I  (Ai),  given  the  values  at  the  four 
integration  points  (Ai,  A2,  A3,  and  A4)  shown  in  Figure  7.3,  are  given  as 

4=14*,  P-D 

1=1 
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Shell 

Stresses 


Thrust  /  foot  =  (<rnrvsl)(  Thickness) 


Figure  7.2.  Calculate  moment  and  thrust  from  shell  stresses 
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where 


and 

r  =  a/3£  ,  5  =  V3/7 


Nx  =0.25(1  -rXl-  5) 
N2  =0.25(1  +  r\l-  s) 
N3=  0.25(1 +  r\l  +  s) 
N4  =0.25(1  -r\l  +  s) 


7.2  Stress  Contour  Calculation 

Stresses  can  also  be  calculated  at  selected  cross  sections  in  the  model.  Cross 
sections  will  pass  through  both  shell  elements  and  superelements.  The  stress 
calculation  is  different  for  each  case.  When  a  section  is  cut  along  a  numbered  or 
lettered  section  line,  the  shells  are  cut  though  the  midplane.  The  shell  element 
stresses  at  the  midplane  are  derived  from  the  shear  and  thrust  values  calculated 
earlier  by  dividing  the  traction  per  unit  length  by  the  thickness.  These  values  are 
transformed  to  the  global  coordinates  using  the  direction  cosines  at  the  node. 
Moment  values  do  not  affect  stresses  at  the  midplane. 

When  sections  are  cut  normal  to  a  shell  element,  stresses  are  also  derived 
from  the  shear,  moment,  and  thrust  values,  but  the  two-dimensional  (2-D)  cross 
section  must  be  constructed  from  the  shell  element  nodal  coordinates  and  the 
thickness  information.  A  pseudo-element  is  constructed  with  the  shell  nodes  on 
the  midplane  and  new  nodes  on  the  top  and  bottom  surfaces  at  each  shell  node,  as 
shown  in  Figure  7.4.  Stress  values  at  these  nodes  are  calculated  by,  essentially, 
inverting  the  procedure  used  earlier  to  calculate  the  shear,  moment,  and  thrust 
from  the  nodal  stresses. 

Where  cross  sections  are  cut  through  superelements,  the  strains  and  stresses 
are  calculated  from  the  nodal  deflections  of  the  original  solid  model  and  the 
strain-displacement  matrices  of  the  hexahedral  and  prism  elements.  The  nodal 
deflections  are  calculated  by  transforming  the  deflections  of  the  shell  nodes 
calculated  in  the  global  finite  element  model,  as  described  above.  The  cross 
sections  are  forced  to  fall  on  element  boundaries,  so  the  strains  and  stresses 
calculated  at  the  nodes  of  the  solid  element  become  nodal  values  for  the  2-D 
quadrilateral  element  that  appears  in  the  cross  section. 


7.3  Shear,  Moment,  and  Thrust  Diagrams 

Once  stresses  on  a  cross  section  have  been  determined,  shear,  moment,  and 
thrust  can  be  calculated  on  any  plane  through  the  section  by  numerically 
integrating  the  stresses.  The  appropriate  stress  values  are  calculated  at  all  element 
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boundaries  intersected  by  the  plane  and  assumed  to  vary  linearly  across  the 
element. 
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Figure  7.4.  Pseudo-element  for  evaluating  shell  stresses 


For  example,  Figure  7.5  shows  the  extensional  stresses  in  the  horizontal 
direction  for  a  typical  section.  Three  finite  elements  were  intersected  by  the 
vertical  line  at  X  =  3.7.  The  thrust  in  the  X  direction  is  determined  by  calculating 
the  net  area  enclosed  by  the  stress  plot.  The  stress  accounted  for  by  the  thrust  is 
then  subtracted  from  the  extensional  stresses.  The  approximate  neutral  axis  can 
be  determined  from  these  adjusted  values.  Moment  is  then  calculated  by  adding 
the  contribution  of  each  area  of  the  stress  plot  multiplied  by  the  distance  to  the 
neutral  axis.  Values  can  be  calculated  at  several  points  along  the  horizontal  axis 
to  construct  a  moment  diagram. 


FE  Model  Stresses 
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Figure  7.5.  Calculation  of  thrust  and  moment  from  cross-section  stresses 
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8  Soil/Pile  Modeling 


8.1  Modeling  Assumptions 

During  initial  fabrication,  the  structure  is  supported  on  the  soil  in  the  dry 
dock.  After  the  structure  is  placed  on  its  underwater  foundation,  it  is  supported 
by  a  combination  of  piles  and/or  drilled  piers,  and  the  grout  pumped  under  the 
structure.  The  nature  of  the  support  condition  will  affect  the  response  of  the 
structure.  The  soil  or  grout  is  modeled  as  an  elastic  foundation  under  the  shell 
elements.  The  supporting  piles  or  drilled  shafts  are  assumed  to  be  rigid  in  the 
vertical  direction,  but  lateral  stiffness  can  be  estimated  and  accounted  for  in  the 
DBA  analysis. 


8.2  Soil  Stiffness 

The  formulation  for  modeling  shells  on  an  elastic  foundation  is  summarized 
in  Chandrupatla  and  Belegundu  (1997).  The  formulation  assumes  that  the  soil 
response  is  linear  elastic.  The  designer  must  know  the  stiffness  of  the  soil  in 
place,  that  is,  the  amount  the  soil  deflects  vertically  when  subjected  to  a  uniform 
pressure.  The  units  for  this  stiffness  are  pressure  divided  by  distance,  for 
example,  100,000  psf/ft  (or  100  kips/ft3).  This  property  is  related  to  the  modulus 
of  vertical  subgrade  reaction,  Kv,  described  by  McCarthy  (1998).  Values  for  a  1- 
sq  ft  plate  range  from  100  to  700  kips/ft3  depending  on  the  soil  type.  The  actual 
settlement  is  a  function  of  the  size  of  the  foundation  and  is  greater  for  larger 
foundations — approaching  4  times  the  settlement  of  the  1-  sq  ft  plate  with  equal 
pressure.  The  effects  of  time-dependent  response  should  also  be  considered  when 
selecting  this  value. 


8.3  Pile  Stiffness 

Methods  for  calculating  pile  stiffness  have  been  developed  by  the  Corps  of 
Engineers’  CASE  Project  (Office,  Chief  of  Engineers  1983).  These  methods  are 
based  on  the  stiffness  of  the  actual  pile  or  drilled  pier,  the  stiffness  of  the 
surrounding  soil  or  bedrock,  and  the  type  of  interaction  assumed  between  the  pile 
and  the  soil. 


Chapter  8  Soil/Pile  Modeling 
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Appendix  A 
Design  Algorithms 


[Reference:  ACI  318-02,  Chapters  13  and  14] 

A.  Flexure  Design  of  Slabs: 

1 .  This  procedure  applies  for  design  of  slab  systems  reinforced  for 
flexure  in  more  than  one  direction,  with  or  without  beams  between 
supports. 

2.  This  procedure  excludes  one-way  slabs  and  soil-supported  slabs. 

3.  Moments  are  determined  using  finite  element  analysis. 

4.  The  maximum  specified  fy  is  60,000  psi  in  shells  and  plates. 

5.  Determine  all  the  data  for  the  analysis:  f ,  Fc,  b,  L. 

(Initially  assume  b  =  1  ft  to  compute  area  of  reinforcement  Asper  foot  of  slab.) 

¥  6.  Determine  the  minimum  slab  thickness,  h: 

Two-way  slab:  Ratio  long/short  span  *2  (b  <  2L). 

Minimum  thickness  for  the  slab  supported  with  beams  on  all 
sides,  no  interior  beams.  For  minimum  thickness,  a 
calculation  is  necessary  to  calculate  the  moment  of  inertia  of 
the  slab  and  beams.  ACI  318  specifies  that  where  beams 
include  portions  of  the  slab  as  flanges,  the  length  of  the 
flange  has  to  equal  hw  or  4*h  (in.),  whichever  is  smaller,  as 
shown  in  the  following  figure: 
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BEAM 


SLAB 


BEAM 


,i*4c,  T'  •  -  •  '  •  ;  , :  '■ 


I  :  l 


Portion  of  the  slab 
included  as  beam 


^  This  portion  is  considered  in  the  moment  of  inertia  calculation. 
Using  the  moment  of  inertia  of  each  section  (slab  and  beam)  in 
both  directions,  a  is  calculated  for  both  directions  as 

Eb*lb 

Of]  2=  — - - ~ 

’  E  *1 

s  s 

S  where  Eb  and  Es  represent  the  modulus  of  elasticity  of  the  beam 
and  of  the  slab,  respectively,  and  Ib  and  Is  represent  the  moment 
of  inertia  of  the  beam  and  of  the  slab,  respectively.  After  ai  and 
a2  are  determined,  am  is  calculated  as 


*  The  minimum  thickness  depends  on  this  value.  If  ctn,  is 
equal  to  or  less  than  0.2,  use 


/*.„,= 


/  *  0.8  + 


200,000 


where  ln  (in.)  is  the  clear  span  of  the  slab,  and  fy  is  in  pounds 
(force)  per  square  inch.  For  am  greater  than  0.2  but  not  greater 
than  2.0,  use 


A.-u  =■ 


/  *  0.8  + 


200,000 


36  +  5.  p{am -0.2) 
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and  for  am  greater  than  2.0,  use 


7.  Cover  and  effective  depth,  d.  Assume  bars  from  following: 

•S  For  spans  up  to  20  ft:  generally  use  No.  4  bars. 

•S  For  spans  up  to  25  ft:  generally  use  No.  5  bars. 

S  For  span  over  25  ft:  generally  use  No.  5-6  bars. 

S  The  clear  cover  (cc)  for  concrete  exposed  to  weather 
conditions: 

•  For  No.  5  bars  and  smaller:  minimum  cc  =  1 .5  in. 

•  For  larger  bars,  minimum  cc  =  2  in. 

■S  Using  these  assumptions,  d  can  be  estimated  for  various  bar 
diameters  (db)  as 

•  For  the  long  span  of  a  flat  plate  or  slab: 

d  =  h  -  cc  -  1.5*db 

•  For  the  short  span  of  a  flat  plate  or  slab: 

d  =  h  -  cc  -  0.5*db 

8.  Compute  a  and  the  reinforcement  area  As  (slab  with  beams): 


Effective 
depth,  d 


T  =  As*  fy 

C  =  0.85*  f  c  *  a  *  b 


(1) 

(2) 


Mn  =  T  (d  -  a/2)  =  As  *  fy  (d  -  a/2)  (3) 
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(4) 


^  Solving  for  a  in  Equations  1  and  2  by  equilibrium  T  =  C: 
a  =  As*fy/(0.85*  fc*b) 

^  Solving  for  As  in  Equation  3: 

As  =  Mn/[fy*  (d  -  a/2)]  (5) 

S  Assuming  an  a,  calculate  A,  from  Equation  5. 

S  With  the  value  of  A,  obtained,  calculate  a  from  Equation  4. 

^  Continue  the  iterations  (calculating  A,  and  a)  until  the 
calculated  a  from  Equation  4  is  almost  equal  to  a  calculated  in 
the  last  iteration. 

9.  Check  As(min): 

^  A,  (mj„),  the  minimum  reinforcement  for  the  slab,  should  be 

at  least  the  minimum  amount  for  shrinkage  and  temperature 
reinforcement: 

•  Slabs  where  Grade  40  or  50  deformed  bar  is 

used . 0.0020 

•  Slab  where  Grade  60  deformed  bars  or  welded  wire 

fabric  (plain  or  deformed)  is  used . 0.0018 

•  Slabs  with  reinforcement  yield  stress  exceeding 
60,000  psi  measured  at  a  yield  strain  of  0.35  %  is 

used . (0.0018  *  60,000)/f 

Note:  the  ratio  of  reinforcement  area  to  gross  concrete  area  cannot  be  less 
than  0.0014. 


•  If  A,  >  As  p  . 

5  s  (min)’  r  A'  mm 

•  If  not,  use  A,  i_._?  as  Af,  and  check  Mn  for  that  As 
(Step  8). 

•  Calculate  a  if  At  changes,  as 

a  =  (A  *  f  )/(0.85  *  r  *  b) 
and  the  new  p : 


A4 
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S  If  As  does  not  change,  calculate 
p  =  As/(b  *  h) 

10.  Check  if  p  is  less  than  0.75  p  bal  (fy  and  f  c  in  psi): 

✓  pbal  =  [(0.85*  ,  *  rc)/f]  *  [87,000/(87,000  +  f )]  \ 

•  for  concrete  strength  f  c<  4000  psi,  ,  =  0.85 

•  for  f  c  between  4,000  and  8,000  psi, 

,  =  1.05  -  0.05*(f /1000) 

•  for  f  c  greater  than  8,000  psi,  ,  =  0.65 
•S  Calculate: 

0.75  p  bal 

S  For  p  <  p  bal ,  f5  =  fy,  and  a  is  sufficient.  If  p  >  p  bal,  the 

slab  would  fail  in  compression  with  fs  >  fy.  Start  the 
procedure  again  and  recalculate  a  (increase  b  or  d). 

S  Check  if  p  is  less  than  0.75  p  bal: 

•  If  p  <  0.75  p  baJ,  p  is  less  than  the  maximum  value 

allowed  by  ACI. 

11.  Compute  Mn  and  Mn: 

✓  Mn  =  A*f  *(d-a/2) 

S  (p  Mn  =  0.9*Mn 

S  Check  with  the  factored  ultimate  moment,  Mu: 

(Mu  was  determined  from  finite  element  analysis) 

•  If  (p  Mn  >  Mu,  proceed. 

•  If  <p  Mn  <  Mu,  stop.  Asis  too  small,  increase  As  and  start 
again. 

12.  Determine  the  As  total  by  multiplying  it  by  b  (total  length  in  feet). 

13.  With  the  total  As  and  the  bar  diameter,  select  the  numbers  of  bars  (n) 
to  use  in  the  design. 

14.  Spacing  for  the  reinforcement  A,: 

S  The  spacing  for  the  reinforcement  A,  for  the  two-way  slab  at 
the  critical  sections  shall  not  exceed  2*h 
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'S  Maximum  spacing  for  shrinkage  and  temperature 
reinforcement  bars: 

•  Shall  not  exceed  3*h  (thickness  of  slab,  inches) 

•  Shall  not  exceed  18  in. 

15.  Calculate  the  mean  spacing  in  the  slab  by 

S  =  b  (total,  in.)/n  (S  <  2*h) 

16.  Repeat  the  procedure  for  maximum  positive  and  negative  moments. 


Negative  moment  Positive 

moment 


17.  For  negative  moment,  d  is  calculated  in  the  same  way  as  for  positive 
moment. 


B.  Shear  Design  of  Slabs: 


1.  Determine  the  shear  from  finite  element  analysis. 

2.  For  shear  design,  it  is  necessary  to  consider  the  behavior  of  slabs  failing 
in  one-way  shear  and  in  two-way  shear. 

3.  It  is  important  to  determine  the  location  of  the  critical  perimeter. 
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4.  One-way  shear: 


(use  the  greater  d  from 
positive  or  negative 
moments  in  each  direction) 


.  Vc  =  2Vfcb„d 

5.  Two-way  shear: 

■/'  Take  the  smaller  of 

•  Vc  =  [2  +  (4/pc)]*rcb0d 

•  Vc  =  (sd)/(b0+2)«fcb0d 

.  Vc  =  4«rcb0 

•  Where 

Pc  =  ratio  of  long  side  to  short  side  of  column, 
reaction  area  or  concentrated  load. 
s  =  30  for  edges  columns 
b0  =  perimeter  of  the  critical  section 

6.  Calculate  the  allowable  shear  for  both  two-  and  one-way  shear: 

(Vu  was  determined  from  finite  element  analysis) 

•  V  =  0.85*V_ 

7.  Check  for  two-way  shear  and  for  one-way  shear: 

•  Check  with  the  factored  ultimate  shear,  Vu: 

•  If  (pVc>  Vu,  proceed.  If  <pVc  <  V„>  stop.  The  shear 
capacity  can  be  increased  by  increasing  the  bo,  and 
by  adding  shear  reinforcement.  Use  closed  stirrups 
with  bars  in  all  four  comers.  The  stirrups  are 
designed  as  for  a  beam,  and  must  be  extended  at 
least  as  far  as  the  critical  section  on  which  cpVc  >  Vu 
(one-way  shear).  These  stirrups  are  not  practical  for 
slabs  less  than  about  12  in.  thick.  Two  different 
sizes  of  stirrups  are  required  to  fit  between  the 
layers  of  flexural  reinforcement  in  the  two 
directions.  They  must  be  spaced  at  d/2,  which 
requires  a  large  number  of  stirrups. 
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C.  Design  of  Walls: 


1.  Minimum  thickness  of  walls: 

^  Exterior  basement  walls  and  foundation  walls:  should  not  be 
less  than  7.5  in. 

^  Bearing  walls  should  not  be  less  than  4  in.  or  L/25,  where  L 
is  the  supported  height  or  length,  whichever  is  shorter. 

2.  Reinforcement  in  walls: 

^  Minimum  vertical  and  horizontal  reinforcement  ratio: 

•  Vertical: 

(a)  Deformed  bars  not  larger  than  No.  5  with 

yield  strength  not  less  than  60,000  psi... 0.0012 


(b)  Other  deformed  bars . 0.00 1 5 

•  Horizontal: 

(a)  Deformed  bars  not  larger  than  No.  5  with  yield 

strength  not  less  than  60,000  psi . 0.0020 

(b)  Other  deformed  bars . 0.0025 


3.  Calculate  the  minimum  area  of  reinforcement  in  each  direction: 

'  A s  min_  vertical  =  ^min  *  k  *  B  for  Vertical  direction 

^  As min_ horizontal  =  ^min  * h* L  horizontal  direction 

4.  Calculate  the  axial  capacity  load  as 

v'  (pPn  =  0.8  *  *  |o.85  *  f  c  *  [a  -A  )+/  *A  \ 

u  '  g  st'  y  st  J 

^  Ag  is  the  gross  area  of  section  in  that  direction,  in.2 
^  Ast  is  the  reinforcement  steel  in  that  direction,  in.2 

5.  Check  if  <p  Pn  >  Pu.  If  not,  increase  Ast  or  h. 

6.  Calculate  the  shear  capacity  of  the  wall  (kips): 

S  d  —  0.8  *  lw  (lw  in  inches).  ACI  uses  this  value  as  the 
effective  depth  for  horizontal  shear  forces. 

^  If  Vc  <  Vu ,  shear  reinforcement  is  needed. 
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7.  Shear  reinforcement  shall  be  provided  as 

A  *  f  *d  i — r 

^  r,--^ — <s*47c*K*d 

si 

y  A  _fas*jA*h*s)^50*bw*s2 

v  f  f 

J  y  J  y 

S  <pV„=<p*(Vc+Vs)<lO*y[frc*h*d 

S  Check  if  V„  >  Vu .  If  not,  add  more  shear  reinforcement. 

8.  Spacing  limit  for  vertical,  horizontal,  and  shear  reinforcement,  s: 

S  Shall  not  exceed:  d/2  <  24  in. 

•/  Horizontal  shear  reinforcement  spacing:  shall  not  exceed 
lw/5,  3*h  or  18  in.  (lw  is  the  horizontal  length  of  wall). 

S  Vertical  shear  reinforcement  spacing:  shall  not  exceed  lw/3, 
3*h  or  18  in. 

9.  Development  length  for  bars  and  hooked  bars  for  each  direction: 

V  The  greater  of  12db  or  d  (effective  depth  at  each  direction). 
^  Hook  length  will  be  12db. 
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